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Abstract. We consider internal shocks as the main dissipation mechanism responsible for the emission in blazars and show 
that it can satisfactorily account for the properties of all blazars. In particular, we extend previous work (Spada et al. 2001) on 
powerful objects, to intermediate (BL Lac) and low power sources (Mkn 421), in order to reproduce the whole of the blazar 
sequence. The model self-consistently treats the dynamics, spectral emission and its variability. The key parameters driving the 
phenomenological sequence are the jet power and the properties of the broad line region, which regulate the cooling efficiency 
of the emitting particles and thus the shape of the spectral energy distribution. By assuming that the remaining parameters 
are similar for all objects it has been possible to reproduce the full range of the observed spectral "states". A more detailed 
comparison of the variability properties shows (for Mkn 421) a good agreement in the X-ray band, while in the optical the 
simulated flux appears to be too variable. For BL Lac lags (~ 10 days) are predicted between the 7-rays and the infrared 
emission. 
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1 . Introduction region. However key issues in the understanding of relativis- 

H tic jets remain open, most notably the jet energetics and the 

. , The discovery that blazars are strong 7-ray emitters together ^^-^^^ acceleration process(es). In order to explore these is- 

with the results of multi wavelength campaigns have allowed ^^^^ ^^^^-^ relationship we have (quantitatively) considered 

to deepen our knowledge on these objects. Their Spectral ^ ^^^^^^-^ ^j^i^j^ ^j^^ pl^^^^^ conditions and their variabil- 

Energy Distribution (SED) is characterized by two broad peaks -^^ ^^^^^^^ f^^^ parameters, but follow from the jet 

(Fossati et al. 1998) strongly variable on different timescales dynamics, thus relating the observed emission properties with 

(Wagner & Witzel 1995; Ulrich, Maraschi & Urry 1997). ^j^^ ^^^^gy transport along the jet. Such a scenario assumes that 

Two main radiation processes dominate the emission, internal shocks are responsible for the dissipation within jets 

namely synchrotron at low frequencies and - plausibly - in- (Rees 1978 and then mostly explored for Gamma-Ray Bursts, 

verse Compton at high energies (see e.g. Sikora 1994 for a re- Rees & Meszaros 1994). 
view). The relative luminosity in the two peaks and their peak 
frequency appear to be functions of the total power (Fossati et 

al 1998, Ghisellini et al. 1998), resulting into a sequence for the The key assumption of the model is that the energy is chan- 
whole of the blazar population properties, ranging from pow- neled into jets in an intermittent way by the central engine, 
erful, low frequency peak, through intermediate, to low power, though such a time dependent process cannot be inferred from 
high frequency peak (blue) blazars (see however Padovani et first principles. Different parts of the jet ('shells') moving at 
al. 2003). The blazar emission is variable on energy dependent different speeds can collide giving rise to shocks and dissipa- 
timescales, typically of weeks-months in the radio and of the tion as non-thermal radiation. The mechanism has a limited 
order of a day in the 7-ray band. efficiency (unless the contrast in Lorentz factors between dif- 
Several studies, mainly based on the modeling of the SED, ferent shells is extremely large, see Beloborodov 2000; Guetta, 
consistently derived the physical parameters of the emitting Spada & Waxman 2001, but also Ghisellini 2002), which has 
to be indeed the case for blazars since most of the energy prop- 
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Beside the low efficiency, the internal shock scenario can 
naturally account for other blazar properties. It predicts that 
jets become radiative at >10^^ — 10^'' cm from the central 
engine, implying a minimum distance from the accretion disk 
and a minimum dimension for the 7-ray source, as required 
by observations (in order to avoid copious pair production; see 
Ghisellini & Madau 1996). Furthermore, successive collisions 
taking place at larger distances have reduced efficiency - since 
the Lorentz factor contrast of the colliding shells decreases 
- explaining why the jet luminosity decreases with distance. 
Finally, the internal shock scenario appears a promising non 
steady state model to account for the observed large amplitude 
variability. 

A detail study of the predictions of this model via numer- 
ical simulations has been carried out by Spada et al. (2001; 
hereafter SOI) for powerful blazars (specifically 3C 279). As 
the SED of this object does not represent the whole of the 
blazar family mentioned above, we focus here on two well 
studied sources representative of extremely blue and interme- 
diate blazars, namely Mkn 421 and BL Lac itself. Rather than 
reproducing in detail particular spectra of specific objects, the 
aim of this work is to determine whether the proposed scenario, 
including the dynamics and emission properties, i) can account 
for the different SED along the blazar sequence, ii) under which 
hypothesis this can occur and iii) whether the internal shock 
model can reproduce the observed variability behavior. 

The outhne of the paper is the following. In §2 we describe 
the hypothesis on the wind dynamics and on the radiative prop- 
erties of the shocked plasma. The results - specifically referring 
to the sources Mkn421 and BL Lac - are presented in §3. In §4 
we draw our conclusions. 

2. Internal shocks and radiated spectra 

2. 1. Physical scenario and dynamics 

In this work we use the approximate model of an unsteady wind 
described in detail in SOI and in the following we report only 
its most relevant assumptions. 

We consider a compact source, of typical dimension Rq ~ 
10^'* cm, which produces an unstable relativistic wind charac- 
terized by an average luminosity Lw The emission from the 
wind is obtained by adding pulses radiated in a series of inter- 
nal shocks that occur in the outflow (Daigne & Mochkovitch 
1998; Spada, Panaitescu & Meszaros 2000; Guetta, Spada & 
Waxman 2001). 

The wind flow is described as a sequence of iV = t^j/U 
shells, where t^f is the overaU duration of the wind emission 
and is the average interval between two consecutive 

shells. Each shell is characterized by four parameters: 

1. the ejection time tj, where the subscript j denotes the j*^ 

sheU, 

2. the Lorentz factor Fj , 

3. the mass Mj, 

4. the width Aj as measured in the lab frame. 

The Lorentz factors of the shells are randomly extracted 
from a uniform distribution between and Tm- The masses 



are also randomly extracted from a uniform distribution with an 
average value Mt/N, where the total mass in the wind Mt is 
normalized to the jet energy budget. The time intervals fj+i —tj 
are drawn randomly from a uniform distribution with an aver- 
age value ty. The initial shell width is of the same order 
of -Ro and the initial internal energy Uj is negligible (all of the 
energy is in the form of bulk kinetic energy). 

The dynamics of the wind expansion is characterized by a 
series of two-shell collisions in which the faster shells catch 
up with the slower ones. The model computes the collision 
time of two shells and replaces them with a new shell having 
mass, Lorentz factor, and internal energy given by the conser- 
vation laws for inelastic collisions (for details see SOI). Before 
colliding, the shells are assumed cold with constant thickness 
A = i?o but after collision they get heated and thus assumed to 
expand at the (comoving) sound velocity. Adiabatic losses dur- 
ing the expansion are taken into account but the bulk Lorentz 
factor is kept constant since the internal energy is always much 
smaller than the bulk kinetic one. The magnetic field is gen- 
erated at each collision without keeping memory of any field 
generated in previous colUsions. 

In each collision a forward (FS) and a reverse (RS) shock 
are formed and propagate into the front and the back shell, re- 
spectively. The plasma parameters behind each shock are de- 
termined by the jump conditions in mildly relativistic shocks, 
namely continuity in the energy density and velocity across the 
contact discontinuity separating the two shells (Panaitescu & 
Meszaros 1999). We also take into account the compression of 
the merged shell width due to the propagation of the FS and 
RS. 

The energy released in each shock is assumed to be dis- 
tributed among protons, electrons and magnetic field with frac- 
tions Cp, Ce and eb. respectively. The resulting magnetic field 
strength allows us to evaluate the synchrotron and the inverse 
Compton emission. 

Electrons are assumed to be accelerated at relativistic 
energies with a power law spectrum, A^(7), above a (ran- 
dom) Lorentz factor 7b. 7b is determined by requiring that 
/ N{'y)^meC^dj equals the fraction Ce of the total available 
internal energy and accounting for the fact that only a frac- 
tion, (e, of electrons is effectively accelerated. We assume that 
below 7b electrons are accelerated with smaller efficiency, cor- 
responding to a flatter energy distribution (with respect to the 
slope above 7b), assumed to be oc 7"^. 

In conclusion, by considering the wind dynamics and shock 
hydrodynamics, the energy density in non-thermal electrons 
and in magnetic field and the typical random Lorentz factor 7b 
can be determined for each collision. 

2.2. Local spectra 

We followed the same basic assumptions of SOI, which we 
briefly recall: 

- The emitting zone is homogeneous with comoving volume 

V — TTtp'^ Rl^ AR' , where ip is the half-opening angle of a 
conical jet and AR' is the comoving thickness of the shell 
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at the collision time, whose expansion while emitting is ne- 
glected for simpUcity. 

The relativistic particles, embedded in a tangled magnetic 
field, have the same energy distribution throughout the re- 
gion corresponding to the shell-shell interaction. This sim- 
plification is justified as we are interested in the "average" 
spectrum, but does not account for spectral details occur- 
ring on a timescale shorter than the light crossing time. 
Even in the case of low power BL Lacs, we consider 
the presence of soft photons external to the jet (identi- 
fied with the emission from the broad line region, BLR). 
The external photon luminosity Lext = aidisk is pro- 
duced within i?BLR and abruptly vanishes beyond, and cor- 
responds to a comoving radiation energy density J7ext = 
(17/12)aLdiskrV(47ri?|LRc) (e.g. GhiselUni & Madau 
1996). 

High and low power blazars have different 
line/photoionizing disk luminosities. Following Kaspi 
et al. (2000), -Rblr and Ldisk are related by -Rblr oc L''^-^^ 
with b ^ 0.7. Thus blazars with weaker broad emission 
lines should have smaller BLR, implying that the first 
collisions occur preferentially outside i?BLR- 
Relativistic electrons are injected with a broken power-law 
energy distribution, oc below 7b and oc 7^* between 
7b and 7niax, with a power L^. The equihbrium particle 
distribution A'' (7), resulting from the continuous injection 
and cooling processes, is determined through the following 
procedure. 

Above the (comoving) random Lorentz factor 7coob cor- 
responding to a cooling time comparable with the shock 
crossing time tcross. the distribution is assumed to have a 
power law index p = s + 1. Below 7cooi there are two dif- 
ferent behaviors: 

i) if 7cooi > 7b, we assume N{'-f) oc 7"** between 7b and 
7cooi, and A^(7) oc 7"^ for 7 < 7b; 

ii) if 7cooi < 7b, ^(7) 0^ 7~^ between 7cooi and 7b, and 
A''(7) oc 7~i for 7 < 7cooi. 

The normalization of A'^(7) is determined according to 
whether particles with 7 ^ 7b can or cannot cool in 
the timescale Across (fast and slow cooUng regime, respec- 
tively): 

Fast cooling regime — If most of Lc is radiated in a 
timescale < Across (i-e- 7b > 7cooi) we apply a luminos- 
ity balance condition: 



Le = VrUeC^ j 7-/V(7)rf7, 



(1) 



where 7 is the radiative cooling rate, including synchrotron, 
synchrotron self-Compton (SSC) and external Compton 
(EC) losses. In order to allow for multiple Compton scatter- 
ings we a priori calculate how many scattering orders nic 
can take place before the Klein-Nishina regime is reached. 
By considering 7b as the relevant energy 



niG 



ln(7b/xB) 
ln(47^/3) 



where xb = hv-Q / {mf.c?') and v-q = eB/{2TrmeC) is 
the (non-relativistic) Larmor frequency. We define Ue = 
Lei {ttR'^c) and introduce the Comptonization parameter y: 



y = ^utAR' j 7^7V(7)d7. 



(3) 



The equality between the injected and the radiated power 
in the fast cooling regime translates into an equality be- 
tween the electron and the radiation energy densities, where 
the latter comprises the synchrotron and the nic inverse 
Compton orders energy density, i.e. 



y + y^+ ....y" 



(4) 



This provides the normahzation of N{'-f), assuming the 
above slopes. 

Slow cooling regime — If most of the power is not radiated 
in fcross> the electrons retain their energy and 



Eo — Lgtn 



I 



c^V / 7V(7)7rf7. 



(5) 



(2) 



Assuming the same distribution for N{'y), its normalization 
can be determined. 

- The slopes of the N{-y) distribution are fixed, but the cool- 
ing energy is found by iteration as 7cooi depends on the 
amount of SSC radiation, which in turn depends on the 
exact shape of N{-y) (because of Klein-Nishina effects). 
Since 7cooi determines the normalization of A?' (7), we re- 
iterate the procedure described above until it converges 
(usually 3-4 iterations are enough). 

- The synchrotron-self-absorption process, the inverse 
Compton emission and the beaming of the radiation are 
treated as in SOI. 

2.3. Observed spectrum 

Both the particle and radiation local spectra thus determined 
are not steady, but roughly correspond to those attained after 
a time Across after the shell-shell interaction, at the maximum 
bolometric radiative output. This spectrum - adopted for the 
entire duration of the emission - has a normalization modulated 
as a function of time. 

Each photon pulse - with its corresponding spectrum - 
starts at the time of a collision and lasts for a duration set by 
the combination of the radiative, shock crossing and angular 
spread timescales. The shape of the pulses (raise and decay in 
the light curves) is assumed symmetric and linear with time 
if its duration is determined by geometrical effects and/or by 
the shock crossing time, as generally occurs for collisions in- 
side the BLR (as typical of powerful blazar like 3C279). If the 
pulse duration is instead set by cooling, the pulse still rises lin- 
early, but decays exponentially [oc exp(— t/3tcooi(7b))], im- 
plying that the rise and decay times can be different. This usu- 
ally occurs in weaker blazars, like Mkn 421. The pulses of in- 
termediate BL Lac objects can be characterized by either of 
these pulse shapes, depending on the broad lines intensity and 
thus on whether colUsions take place within .Rblr- 
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Fig. 1. The SED of Mkn421 during different observational 
campaigns, as labeled. For the data points see Costamante & 
Ghisellini (2002) and references therein. 



o 

M 



O 
,-J 



-10 - 



■11 



■12 



■13 « 



BL Lac 




- Q ■ 1999 June O 1995 Jan-Feb 

o 



1997 Flare X Not sim. data 



O 1995 Nov 



10 



15 20 
Log V [Hz] 



25 



Fig. 2. The SED of BL Lac during different observational cam- 
paigns, as labeled. For the data points see Ravasio et al. (2002) 
and references therein. 



The observed spectrum is determined by all of the photons 
reaching the observer simultaneously: since these are produced 
at different distances, their different light propagations times 
should be taken into account. In other words, the resulting spec- 
trum is obtained as the convolution of the emission produced 
by shocks simultaneously active in the observer frame, taking 
into account the radius at which collisions occur and the photon 
travel path. 



The main aim of this work is to explore whether the internal 
shock scenario can account for the difference of the SED along 
the blazar sequence. We thus a) concentrate on the global prop- 
erties, rather than to explain details of the spectrum of a spe- 
cific object; b) apply the model to Mkn 421 and BL Lac itself 
which, besides being among the best studied BL Lac objects, 
are representative of extremely blue and intermediate blazars, 
respectively, complementing the study already performed for 
the powerful blazar 3C 279 (SOI). 

Mkn 421 — This is one example of bright, extremely blue 
and strong TeV emitting blazar The peaks of its synchrotron 
and inverse Compton spectra are at ~ keV and a few hundreds 
GeV energies, respectively. The bolometric (observed) power 
is of order 8 x IC*^ erg s^^. The SED of coordinated simul- 
taneous observational campaigns of Mkn 421 are reported in 
Fig. 1 (Costamante & Ghiselhni 2002). 

BL Lac — The prototype of the BL Lacertae class is an 
intermediate blazars (e.g. Fossati et al. 1998). It is a 1 Jy source, 
with broad band peaks in the near IR and MeV-GeV band and 
its total bolometric power is of order ^ 10^^ erg s~^. Despite of 
BL Lac being the prototype of the BL Lac class, it sporadically 
showed optical emission lines (with maximum observed EW 
- 6 A; Vermeulen et al. 1995; Corbett et al. 2000). In July 1997 
it underwent a major outburst, followed in the optical and by 
the CGRO (EGRET), RossiXJE and ASCA satellites. During 
this flare the entire SED dramatically changed, as illustrated in 
Fig. 2, with EGRET detecting a flux (>100 MeV) a factor 3.5 
times higher than that observed in 1995 (Bloom et al. 1997). 

In order to reproduce the SED and variability of these two 
sources, we took as an initial set of physical parameters that 
determined for 3C 279 (see SOI), and changed as few as pos- 
sible of them to account for the different broad band SED and 
'spectral states' of these sources. 

It turned out that the global properties of the blazar "se- 
quence" can be reproduced by only changing two key quanti- 
ties, namely the jet luminosity Ljct and the luminosity Lblr 
and extension i?BLR of the external photon field. The other 
parameters [eb, Ce, the slope n of the (high energy) electron 
distribution, the observing angle 6] are only constrained to fit 
the specific SED. The remaining quantities are at most slightly 
varied just to optimize the fit. In particular eo has been fixed for 
all sources to 0.5. In Table 1 we list all of the input parameters 
adopted to reproduce the time dependent behavior of the simu- 
lated sources (for convenience of the reader we also report the 
same parameters for 3C 279 from SOI). 

Let us now examine the physical quantities inferred from 
the modelling and their dependence on R. 

Efficiency — The radiative efficiency represents the fraction of 
the dissipated energy at each location channeled into relativistic 
electrons. In absence of non-radiative losses this energy should 
be all radiated. As shown in Fig. 3 and 4 for Mkn 421 and BL 
Lac, respectively, the efficiency decreases simply due to the de- 
crease in the contrast of bulk Lorentz factors of colliding shells. 
The apparent different efficiencies of BL Lac and Mkn 421 re- 
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Fig. 3. Radiative efficiency as a function of the collision radius 
for Mkn 421. The soUd line refers to the global efficiency, i.e. 
the fraction of the total wind kinetic energy (E^ = Lw x iw) 
radiated on scales smaller than any given radius; the shaded 
histogram shows the differential efficiency, namely the fraction 
of Ey, radiated for a given radius interval. The cone-like insert 
in the upper part of the figure shows a grey-tone representation 
of the differential efficiency of the jet: the darker the color the 
higher the efficiency. 
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Fig. 4. Radiative efficiency versus the colhsion radius for BL 
Lac. The soUd line refers to the global efficiency, i.e. the frac- 
tion of the total wind kinetic energy (E^, = L„ x t^,) radi- 
ated on scales smaller than a given radius; the shaded histogram 
shows instead the differential efficiency, namely the fraction of 
E„ radiated for a given radius interval. The cone-like insert in 
the upper part of the figure shows a grey-tone representation 
of the differential efficiency of the jet: the darker the color the 
higher the efficiency. The vertical Une indicates .Rblr- 



fleet only the different viewing angle (chosen to reproduce the 
detailed properties of the two sources). 

Parameters profiles: T, B, 7b — Fig. 5 and Fig. 6 (top three 
panels) show how the bulk Lorentz factor, the magnetic field, 
the Lorentz factor 7b of the injected electrons evolve with dis- 
tance from the nucleus, for Mkn 421 and BL Lac, respec- 
tively. Different symbols correspond to shells having experi- 
enced a different number of collisions. The bulk Lorentz fac- 
tors tend toward average values as result of successive colh- 
sions and merging of the shells. The range and the asymptotic 
bulk Lorentz factor V are fully consistent with the typical val- 
ues inferred from observations. The magnetic field decreases 
roughly as R~^l'^ for both sources: it is worth remarking that 
we assume that a constant fraction (ee) of the available energy 
is converted into magnetic field and no field component sur- 
vives between successive collisions. The minimum energy of 
the injected electrons, ^fbm,pC?, is also a decreasing function of 
R, due to the decrease in the bulk Lorentz factors contrast, im- 
plying decreasing efficiency and mean energy per particle. 

In general the only difference between Mkn 421 and BL 
Lac is the smaller spread of the parameters for the former, due 
to the fact that in Mkn421 all the coUisions occur in the same 
region. 

Parameters profiles: 7pcak, ^'peak — The two bottom panels of 
Figs. 5, 6 show the spatial profile of the energy of the elec- 
trons emitting at the peak of the SED, mcC^7pcak, and the cor- 
responding observed synchrotron peak frequency, z^peak- 

A significantly different behavior characterizes the two 
sources. In Mkn 421 7peak is a continuous decreasing func- 
tion of R and attains large values. This follows from the fact 
that in low power blazars only high energy electrons can ra- 
diatively cool in a finite injection time, thus determining high 
values of 7peak- Specifically, its value is determined by 7max 
at each collision radius, as set by the spectral fit. The behavior 
of i^peak simply follows that of 7peak and B. On the contrary 
for BL Lac 7peak (and consequently fpeak) shows a disconti- 
nuity. This occurs at /Jblr: inside the BLR electrons are in 
the fast coohng regime and thus 7peak = 76 while outside the 
BLR the behavior is similar to that of Mkn 421, namely the ra- 
diative cooling time exceeds the shock crossing time, implying 
7poak = 7max- Notc that for 3C 279 the BLR is even larger and 
7poak = 76 for almost all of the collisions. 

Spectra at various R — In order to understand how and where 
the overall spectrum is generated, the spectra produced at dif- 
ferent radii are shown in Fig. 7 and Fig. 8: in the top panel 
the average "instantaneous" spectrum of a single shell at R 
is shown while in the bottom panel the spectra are integrated 
for all of the shells emitting at a certain time (the flux is thus 
weighted for the duration of the emission in each shell at a 
given frequency). While the single shell emission is qualita- 
tively similar for the two sources, the integrated spectra of BL 
Lac do not simply reproduce the single shell behavior. This is 
due to the different (shorter) duration of each single shell emis- 
sion, due to the fast cooling of the electrons (caused by the 
external Compton). The end result is that in Mkn 421 most of 
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Fig. 5. Evolution of the post-shock shell parameters as a func- 
tion of the collision radius for Mkn 42 1 . The upper panel shows 
the bulk Lorentz factor of the merged shell after each collision; 
the second panel represents the value of the generated magnetic 
field (eb = 6 X 10^'*). The third panel shows the minimum ran- 
dom Lorentz factor of the injected electrons 7b (for Cc = 0.5; 

= 0.08). In the fourth and last panels 7poak and are 
plotted. Only one tenth of the points is displayed for clarity. 
Circles refer to collisions in which none of the shells has ever 
collided before, stars to collisions in which one of the shells 
has collided, while triangles indicate collisions between shells 
that have both collided before. 



the emission is produced in the inner jet, while in BL Lac the 
inner part only dominates at the highest frequencies of the two 
spectral peaks, with interesting consequences on the correla- 
tions between variability in different bands. 

It is also worth commenting on the extreme differences in 
the SED of BL Lac at different epochs (as shown in Fig. 2). 
This has been already interpreted as due to variable broad emis- 
sion lines in the spectrum (mainly Ha and Hp), which would 
cause a variable contribution of seed photons for the exter- 
nal Compton process. The maximum observed total luminosity 
from the BLR is 5-8 x 10^^ erg s~^ (see e.g. Vermulen et al., for 
Hq = 65 km s~^ Mpc~^). If one adopts the correlation (Kaspi 
et al. 2000) between i?BLR and the ionizing continuum lumi- 
nosity (here assumed to be ~ 10 times Lblr), ^blr ^ 5-7 
xlO^^ cm. This inferred distance is indeed similar to the typ- 
ical distance of the first shell-shell collisions, suggesting that 
the collisions sometimes occurs within i?BLR, sometimes out- 
side it (see also Ravasio et al. 2002, 2003). Dramatic spectral 
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Fig. 6. Evolution of the post-shock shell parameters as a func- 
tion of the collision radius for BL Lac. The upper panel shows 
the bulk Lorentz factor of the merged shell after each collision; 
the second panel represents the value of the generated magnetic 
field (eb = 10^^). The third panel shows the minimum ran- 
dom Lorentz factor of the injected electrons 7b (for eo = 0.5; 

= 10^^). In the fourth and last panels 7pcak and i^°|?ak are 
plotted. Only one tenth of the points is displayed for clarity. 
Circles refer to collisions in which none of the shells has ever 
collided before, stars to collisions in which one of the shells 
has collided, while triangles indicate collisions between shells 
that have both collided before. 



changes, such as those between the 1995 and the 1997 SED 
(see Fig. 2) can be therefore accounted for in this scenario. 

Light curves — The simulations performed aimed at repro- 
ducing (and successfully did) all of the flux range of the 
SED observed in the various multiwavelength campaigns. 
However the model also predicts the detailed variability be- 
havior at all frequencies (see also Tanihata et al. 2003) which 
has not been considered a priori to constrain the model. 
Let us then examine the light curves in different bands, re- 
ported in Fig. 9 and Fig. 10 for the two sources, and sim- 
ply compare their variability properties with observations. 
A movie showing the spectral variability during the whole 
simulation for Mkn 421 (and 3C 279) can be found at the 
URL http://ares.merate.mi.astro.it/'-~^gabriele/421/index.html 
The long term 'periodicity' that appears in the light curves is 
due to the sinusoidal modulation chosen for the mass distribu- 
tions of the shells. 
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Simulated lightcurves 
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Fig. 7. Average spectra of collisions for Mkn 42 1 . The upper 
panel shows the average spectrum of the shells without taking 
into account the duration of the emission. In the lower panel, in- 
stead, the duration of the emission of each shell has been taken 
into account. 
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Fig. 8. Average spectra of collisions for BL Lac. The upper 
panel shows the average spectrum of the shells without taking 
into account the duration of the emission. In the lower panel, in- 
stead, the duration of the emission of each shell has been taken 
into account. 
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Fig. 9. Simulated light curves for Mkn 421 in different bands, 
as labeled. 
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In order to assess the potential of the model to reproduce 
the data, we focused on comparing the simulated light curves 
with the observed ones, by calculating basic indicators, namely 
the r.m.s., the cumulative distribution of the variance [(Flux - 
( Flux ))^], and the minimum variability doubling timescales, 
defined as ( Flux )At/{A Flux). In particular we focused on 
the X-ray and optical bands for Mkn 421. For the optical vari- 
abiUty we sampled the simulated curves with a frequency of 
about once a day for all the 300 days covered by the simulation 
and compared the results with the photometric light curves in 
the F band presented by Tosti et al. (1998). For the X-rays, we 
considered the results presented by Tanihata et al. (2001) on 
6 day continuous observations of the source with ASCA and 
biimed the observed 0.2-2 keV Ught curve over the same time 
bins (0.125 days) of the simulations. 

In the optical band the simulated Ught curves show signif- 
icantly larger variability than the observed one (as from Tosti 
et al. 1998): the r.m.s. from the observations is about 4 times 
smaller than that derived from the model, and the Kolmogorov- 
Smimoff test (KS) on the variance distributions rejects the null 
hypothesis. On the contrary, the comparison reveals promising 
results for the variability in the X-ray band. Here, the observed 
and simulated r.m.s. are comparable (0.2 vs 0.3) and the KS 
test shows that the variance distributions are marginally com- 
patible (5 per cent level). The X-ray variability properties of 
the simulated and observed light curves appear to be similar 
also in terms of the variabiUty timescale, which is of the order 
of 0.3-0.4 day. 

While a proper comparison of the variability properties 
would require a more sophisticated analysis and better sam- 
pled data, these tentative results show that the variability in the 
X-rays is reasonably reproduced by the model. A significant 
discrepancy appears instead in the optical, where the model 
predicts larger variabiUty: this could be in principle improved 



by modifying the model parameters (e.g. reducing the Lorentz 
factor contrast and modifying its distribution), although beyond 
the aim of this paper. We however notice that in the case of Mkn 
421 the possible role of the host galaxy emission in reducing 
the optical variabiUty ampUtude should be assessed. 

Finally, we computed the cross-correlations between the 
variability in different energy bands. The results are presented 
in Fig. 11 and 12 for Mkn 421 and BL Lac, respectively. As 
discussed above for Mkn 421 7pcak is determined, at all radii, 
by 7max, accounting for the absence of lags between the differ- 
ent frequencies (Fig. 11). On the contrary, in BL Lac lags (~ 
10 days) are present between the infrared and the 7-ray band. 
In fact, in this source the 7-ray emission is dominated by the 
external Compton within i?BLR , while the infrared one by syn- 
chrotron radiation from collisions occurring outside the broad 
line region. We remind that in the case of 3C 279 (SOI) the in- 
frared emission is also determined by colUsions outside i?BLR, 
resulting in lags between the IR and the 7-rays of ~ 40 days. 

4. Discussion 

In this work we have considered internal shocks as the dissi- 
pation mechanism responsible for the emission in blazars. As 
mentioned in the Introduction, this scenario is currently the 
most accredited to explain the prompt gamma-ray emission in 
Gamma-Ray Bursts, despite of an efficiency problem (Lazzati, 
GhiselUni & Celotti 1999). In blazars instead the relatively low 
efficiency is in fact required as, at least in powerful objects, 
most of the jet power is carried to the large scales. 

In SOI it was found that the internal shock model was suc- 
cessful in reproducing the observed SED and variability prop- 
erties of a powerful blazar, namely 3C 279. However, the SED 
of this object does not represent the whole of the blazar class. 
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Fig. 11. Results of the cross correlations analysis for Mkn 421, 
between the 7-ray and X-ray (solid grey Une), optical (soUd 
dark line) and infrared (dashed line) light curves. There are no 
delays between the different frequencies. 
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Fig. 12. Results of the cross correlations analysis for BL Lac, 
between the 7-ray and X-ray (solid grey Une), optical (solid 
dark line) and infrared (dashed line) Ught curves. There is a 
delay between the 7-ray and the infrared emission of roughly 
10 days. 
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Table 1. Model input parameters for Mkn 421, BL Lac and 3C279 



which covers a wide range of spectral characteristics (luminos- 
ity, frequency of the peaks of the emission and their relative in- 
tensity). Interestingly these parameters appear to be correlated 
and the whole class can be seen as a sequence: the frequency 
and the ratio of the low vs high energy peak intensity increase 
with decreasing luminosity (Fossati et al. 1998). 

We show here that the internal shock model can satisfacto- 
rily account also for the properties of the lower power blazars. 
The key parameters driving the phenomenological sequence 
are the jet power (proportional to the radiated one) and the in- 
tensity of the broad lines. These parameters in turn regulate the 
SED shape, as they control the cooling efficiency of the emit- 
ting particles. The global radiative efficiency appears instead to 
be similar for all of the sources examined. 

The internal shock scenario determines a characteristic 
time interval for the injection of relativistic electrons of the or- 
der of the dynamical timescale, when the intensity of the spec- 
trum from each collision is maximized. Two regimes are rele- 
vant: fast and slow cooUng, corresponding to whether electrons 
of energy 7b can or cannot radiatively cool in the dynamical 
time. In highly powerful blazars the fast cooling regime dom- 
inates in the inner regions (within the BLR) where also most 
of the power is dissipated. Consequently the peak frequencies 
are produced by electrons of energy 7bmeC^. In the weakest 
blazars instead the slow cooling regime prevails over the whole 
jet. Consequently only the highest energy electrons can cool 
over such timescale: the peak frequencies thus shift to high val- 
ues. 

Between these two extremes there are intermediate sources, 
like BL Lac, with broad lines of intermediate intensity, pro- 
duced at a distance within which a few shell-shell collisions 
can occasionally take place. Observationally this corresponds 
to SED with moderate Compton to synchrotron luminosity ra- 
tio, as the scattered seed photons are only the synchrotron ones 
(in collisions outside the BLR). The rare colUsions within the 
BLR give rise to dramatic changes in the SED characterized by 
a large increase of the Compton component (e.g. BL Lac itself, 
see Fig. 2). 

The model considered here, which considers self- 
consistently the dynamics and spectral emission, predicts also 
the time dependent spectral properties. In general, the selected 
parameters allow to reproduce the full range of spectral 'states' 
observed in both BL Lac and Mkn 421 . A more detailed com- 
parison performed for Mkn 421, and based on the r.m.s. and 
variabihty timescales, shows good agreement for the X-ray 
variabiUty properties, while the simulated optical variations ap- 
pear to be too large with respect to the considered observed 
fight curve. 



An analysis of the predicted cross-correlated variabifity be- 
tween the 7-ray and other bands reveals that only for BL Lac 
lags (~ 10 days) are expected between the 7-ray s and the in- 
frared emission. 

We conclude that the internal shock scenario can account 
for the the main properties of blazars and the 'blazar sequence' . 

While it has been previously pointed out that the key quan- 
tity in reproducing the different characteristics of blazars is the 
ratio between the jet and the disc (i.e. broad lines) luminosi- 
ties, the internal shock scenario discussed here provides more 
physical insights, namely i) directly connects the radiated jet 
luminosity with the jet effective power and ii) accounts for the 
'preferred' distance where most of the luminosity is dissipated. 
While the latter is similar for high and low power blazars, the 
BLR is instead located at different distances in the different 
sub-classes of objects (as determined by the ionizing luminos- 
ity). 

From a more theoretical point of view, the model does not 
(yet) address the role of a seed magnetic field amplified by the 
shell-shell collisions. This would probably lead to a faster syn- 
chrotron cooling on the large scales and a flatter dependence of 
the B field from R (Fig. 3). We intend to further explore this 
issue. Furthermore, only internal shocks have been considered 
in this scenario. It is conceivable that at large jet scales some 
entrairraient may occur causing an interaction of the jet with 
the external medium, possibly leading to external shocks. As a 
result the radiative efficiency and large scales synchrotron and 
inverse Compton emission could be enhanced. Clearly, the in- 
verse Compton emission could be also enhanced if a significant 
photons field is present externally to the jet even on large scales 
(such as microwave background and/or beamed nuclear radia- 
tion and/or dust emission, e.g. Celotti, GhiselUni & Chiaberge 
2001;Sikoraetal. 2002). 

Observationally the improved detection sensitivity of the 
planned TeV Cherenkov telescopes - such as VERITAS, 
HESS, MAGIC - will allow in the near future to measure emis- 
sion from BL Lacs in low luminosity states, and thus to esti- 
mate their flare activity duty cycles. Since in the internal shock 
scenario the TeV emission is largely produced by the few pow- 
erful internal collisions, it will be then possible to further con- 
strain model parameters such as the range of bulk Lorentz fac- 
tors and the initial separation of the shells. Analogous results 
will be likely achieved for the more powerful blazars in the 
GeV band flianks to flie AGILE and GLAST sateUites. 
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